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Abstract 

A search for quark compositeness in the form of quark contact interactions, based 
on hadronic jet pairs (dijets) produced in proton-proton collisions at a/ s = 7 TeV, 
is described. The data sample of the study corresponds to an integrated luminosity 
of 2.9 pb _1 collected with the CMS detector at the LHC. The dijet centrality ratio, 
which quantifies the angular distribution of the dijets, is measured as a function of 
the invariant mass of the dijet system and is found to agree with the predictions of 
the Standard Model. A statistical analysis of the data provides a lower limit on the 
energy scale of quark contact interactions. The sensitivity of the analysis is such that 
the expected limit is 2.9 TeV; because the observed value of the centrality ratio at 
high invariant mass is below the expectation, the observed limit is 4.0 TeV at the 95% 
confidence level. 
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In the Standard Model (SM), most high energy proton-proton collisions are described by the 
scattering of partons (quarks or gluons) in the framework of Quantum Chromodynamics (QCD). 
The outgoing partons manifest themselves as two or more jets of hadrons, with the pseudora- 
pidity tj of the jets depending on the parton scattering angle. In QCD, the jet production rate 
peaks at large \n\ because the scattering is dominated by f -channel processes. Several new 
physics scenarios, including models of quark compositeness, produce a more isotropic angular 
distribution leading to enhanced jet production at smaller values of \n\ QHH- Other models of 
new physics predict the opposite: a decrease in jet production at small \t]\ compared with the 
SM 0. 

The dijet system consists of the two jets with the highest transverse momenta pj in an event 
(the leading jets) with invariant mass m/y. The dijet centrality ratio R )? is defined as the number 
of events with the two leading jets in the region < 0.7 (inner events) divided by the number 
of events with the two jets in the region 0.7 < \n\ < 1.3 (outer events). Because many systematic 
effects cancel in this ratio, R, 7 provides an accurate test of QCD and is sensitive to new physics. 

In this Letter we report a measurement of R )? in proton-proton collisions at y/s = 7 TeV. The 
analysis is based on a data sample corresponding to an integrated luminosity of 2.9 ± 0.3 pb _1 
collected with the Compact Muon Solenoid (CMS) detector at the CERN Large Hadron Col- 
lider. 

The dijet centrality ratio is measured as a function of m;; and is compared with predictions 
from perturbative QCD calculations performed at next-to-leading order (NLO) accuracy with 
the NLOJET++ program HHU in the FASTNLO framework [10[. We also compare the measured 
R,7 with a QCD prediction obtained from the PYTHIA 6.420 event generator [llj with the D6T 
set of parameters [12J. We use CTEQ6.6 parton distribution functions (PDFs) [13j for the NLO 
calculation and CTEQ6LL (HI for the PYTHIA6 prediction. The effect of the CMS detector 
simulation |[T5| on the predictions for R^ is negligible. 

We use Ryf to search for evidence that quarks are composite particles. Quark compositeness 
at an energy scale A would appear at lower energies as a contact interaction, yielding an n 
distribution different from that predicted by QCD. We consider a model of contact interactions 
between left-handed quarks in the process qq —¥ qq described by the effective Lagrangian L qq = 
(±27r/ A 2 )(^ l 7^£)(^ l 7^) d). We choose the positive sign of L™ because it yields a more 
conservative exclusion limit on A (by about 5%) than the negative sign. In QCD, R^ is nearly 
independent of mu, with a value near 0.5. The presence of quark contact interactions described 
by L qq would cause R^ to increase rapidly above a value of m« that depends on A. Previous 
searches for quark compositeness described by the interaction L qq exclude A < 3.4 TeV at the 
95% confidence level (CL) flUHZll- 

A detailed description of the CMS detector can be found elsewhere §22\ . The CMS coordinate 
system has the origin at the center of the detector, the z-axis along the direction of the counter- 
clockwise beam, and the transverse plane perpendicular to the beam axis; (p is the azimuthal 
angle, 6 the polar angle, and n = — In (tan [0/2]) the pseudorapidity. The central feature of 
the CMS apparatus is a superconducting solenoid that surrounds the silicon pixel and strip 
tracker as well as the barrel and endcap calorimeters (covering the region \r/\ < 3): a lead 
tungstate crystal electromagnetic calorimeter (ECAL) and a brass-scintillator hadronic calori- 
meter (HCAL). The ECAL barrel extends to \rj\ = 1.479 and the HCAL barrel to \n\ = 1.305. 
The HCAL and ECAL cells are grouped into towers projecting radially outward from the ori- 
gin. In the region \n\ < 1.74 these calorimeter towers have width An = A(p = 0.087. ECAL and 
HCAL cell energies above noise suppression thresholds are summed within each projective 
tower to define the calorimeter tower energy. 
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We reconstruct jets by applying the anti-Air clustering algorithm [23J to the calorimeter towers 
with the distance parameter R = 0.7. The energy E and momentum p of a jet are defined as the 
scalar and vector sums, respectively of the calorimeter cell energies associated with the jet. We 
apply pj- and ^-dependent scales to E and p to correct for the non-linearity and non-uniformity 
of the calorimeter response. The jet energy corrections and resolutions are determined and 
validated using simulated, test beam, and collision data |24] . 

A set of independent single-jet triggers with varying thresholds on uncorrected jet pj is em- 
ployed in the online trigger system. We use data from three of these triggers, with thresholds 
of 15, 30, and 50 GeV, in the rrijj ranges where the triggers have efficiency greater than 99.5% 
for both inner and outer events. By studying the relative efficiency of parallel triggers in the 
collision data, we determine that these three triggers meet this efficiency requirement for mjj 
greater than 156, 244, and 354 GeV, respectively, where these values are three of the predefined 
bin edges for mjj. The requirement of mjj > 156 GeV results in a minimum jet pj of 25 GeV. 

To remove potential instrumental and non-collision backgrounds we impose the following re- 
quirements: events must have a primary vertex reconstructed with |z| < 24 cm [25]; jets must 
have at least 1% of their total energy detected in the ECAL, no more than 98% of their energy 
detected in a single HCAL photodetector, and no more than 90% of their energy in a single 
calorimeter cell (ECAL or HCAL). These jet identification criteria remove less than 0.1% of the 
events that pass the mjj threshold and rj constraints. 

In Fig. [l] we show the observed numbers of inner and outer dijet events and R^ in bins of mjj) 
the bin widths roughly correspond to the mjj resolution. The event counts, which are corrected 
for the trigger reduction factors (prescales), fall steeply with increasing mjj. We compare R^ 
with NLO and PYTHIA6 predictions for mjj values up to 1120 GeV. The error bars represent the 
combination of statistical and experimental systematic uncertainties (described in detail later). 
The horizontal lines near the end of the error bars denote the statistical uncertainty on this ratio 
of Poisson-distributed variables computed with Clopper-Pearson intervals |26j. 

We apply an m^-dependent correction to the NLO prediction to account for non-perturbative 
effects of hadronization and multiple parton interactions. This correction, which is approxi- 
mately 10% at low ntjj and 2% for mjj greater than 400 GeV, is obtained from PYTHIA6. The 
predictions of PYTHIA6 and HERWIG++ [27| for this correction agree to within a few percent. 

The NLO prediction is shown as a band that accounts for uncertainties related to the choices of 
the renormalization scale p.R, the factorization scale p?, and the PDFs used in the calculation. 
The scale uncertainties, which are approximately 3^% depending on mjj, are evaluated by 
varying the scales from the default choice of /yr = p? = Pt to pj/2, pj, and 2pj in the following 
six combinations: (p R ,p F ) = (p T /2,p T /2), (2p T ,2p T ), (p T ,p T /2), (p T ,2p T ), (p T /2, p T ), and 
{2pj, pj). The PDF uncertainties are estimated with repeated evaluations of the NLO-predicted 
R v for the PDFs in the CTEQ6.6, MSTW2008 |28l, and NNPDF2.0 EU sets and are found to 
be less than 1%. The band also includes the uncertainty arising from the correction for non- 
perturbative effects, which we conservatively take to be 20% of the correction factor. 

The measured R, ? is nearly flat with a value of about 0.5 as predicted by both the corrected 
NLO calculation and PYTHIA6. The observed average ratio is about 7% lower than that of 
the corrected NLO prediction, and about 7% higher than that of the PYTHIA6 prediction. The 
data are in better agreement with the corrected NLO prediction at low mjj, where the signifi- 
cant non-perturbative corrections improve the agreement, and with the PYTHIA6 prediction at 
intermediate and high mu. 

To test for the presence of quark compositeness with R n , we employ a log-likelihood-ratio 
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Figure 1: (a) Event counts corrected for trigger prescales for inner (solid circles) and outer 
(open boxes) dijets and (b) the observed as functions of m«. We compare R n with pre- 
dictions for QCD from PYTHIA6 (dashed line), NLO calculations (dotted line), and NLO plus 
non-perturbative corrections (solid line) and its uncertainty (band). 



statistic that compares the likelihood of the null (QCD only) hypothesis £qcd with that of 
the alternative hypothesis that quark contact interactions are present in addition to QCD C a \ t : 

Tl LL = In £ alt - In £q C d- (1) 

The total likelihood is the product of the individual bin likelihoods, which for ntu bin i is 

£/ = 'P(n to t,i\Htot,i)S(n in/i \n tot j,pi), (2) 

where the first factor is the Poisson probability to observe n tot/ j events when expecting ^ tot ; - and 
the second is the binomial probability to observe n m ^ inner events given n tot/! and a predicted 
probability to be inner of p t ■ (p = R v / (1 + ) ) . Since the first factor in Eq. j2| contains no infor- 
mation on Rjj, we remove it from the statistical inference by conditioning the probabilities by 
the observed values of n tot/1 - ||30[|3T1 . We compare the value of TZn in the data with distributions 
of the expected values for both hypotheses, obtained from ensembles of pseudoexperiments, to 
either claim the discovery of quark compositeness or set exclusion bounds on the composite- 
ness scale A with the frequentist-inspired CL S method |32|. This method provides protection 
against an exclusion claim when the data have little sensitivity to the new physics. 

We use the NLO prediction corrected for non-perturbative effects to describe the shape of R^ for 
the null hypothesis. To minimize the effect of potential discrepancies between the NLO predic- 
tion and actual QCD dijet production, we include an overall offset of R, 7 in the null hypothesis. 
This offset is determined with the data in the nijj range between 490 and 790 GeV. (The lower 
bound is chosen to avoid the region where non-perturbative corrections are significant, and the 
upper bound is chosen to avoid the signal region for compositeness.) As noted above, the data 
lie below the NLO prediction, yielding an offset of &R V = -0.050 ± 0.021 (stat.) ± 0.039(syst.). 
Using ensembles of simulated data, we determine that the probability (p-value) for observing 
| AK^ | > 0.050, given the NLO prediction, is 0.29. 
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PYTHIA6 is used to describe R^ for the alternative hypothesis. We apply an m^-dependent cor- 
rection that accounts for NLO contributions to the QCD part of this prediction. We do not apply 
this correction, which is derived for f -channel QCD processes, to the contact interaction part of 
the prediction because it is not physically motivated and yields less conservative exclusion lim- 
its on A. Since the contact interaction model is not valid for m» near the compositeness scale, 
we exclude data above a A-dependent mu threshold for the testing of each A value hypothesis. 
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Figure 2: The observed dijet centrality ratio as a function of ttiji compared with the null (QCD) 
hypothesis (solid line), including the total systematic uncertainty (band), and to hypotheses of 
quark contact interactions with A = 3 TeV (dotted line) and 4 TeV (dashed line). 



In Table [T] we report the systematic uncertainties related to the measurement of K, 7 and the 
NLO QCD model. The dominant source of uncertainty on the measurement is the 1% uncer- 
tainty in the relative jet energy scale (JES) between the inner and outer rj regions, which results 
in a 5-13% uncertainty on depending on rrijj. This relative uncertainty has a much larger 
impact than a 10% uncertainty on the JES common to both regions. For the QCD model, the 
sources of uncertainty include the choice of scale and PDFs in the NLO calculation and the 
non-perturbative corrections described above. In addition, we take the statistical uncertainty 
on the offset described above and the difference between the PYTHIA6 and NLO predictions 
as systematic uncertainties related to our choice of model. For the compositeness hypothe- 
sis, R^ increases steeply with mu, and the 10% uncertainty on the absolute JES dominates the 
uncertainty on the A scale being probed. 

Figure [2] shows our data in comparison with the null hypothesis. Alternative hypotheses with 
contact interaction scales of A = 3 and 4 TeV are also shown. In this figure, the data from the 
15 sparsely populated mu bins in the range 1530-3020 GeV are combined into a single bin for 
presentation purposes. The band indicates the total systematic uncertainty, which is included 
in the ensembles of pseudoexperiments with the method of Ref. |33|; i.e., the uncertainties 
enter the ensembles as nuisance parameters that affect the expected numbers of inner and outer 
events. 

To quantify the agreement of the data with the SM expectation, we determine the offset of 
the data with respect to the NLO model for the full mu range, finding —0.037 ± 0.007(stat.) ± 
0.039 (syst.) with a p-value of 0.34. Given this consistency of the data with the QCD hypothesis, 
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Table 1: Systematic uncertainties on R^ related to the measurement of R )? (detector uncertain- 
ties) and to the QCD model (model uncertainties). For each source of uncertainty, we show the 
range of values over the entire mu range and at a representative point in the signal region. 



Source 


|— i ii n 

Full Range 


m.jj = 1.6 TeV 


Detector uncertainty 






Relative JES 


0.02-0.05 


0.032 


Absolute TES 


00-0 03 

\J •\J\J \_/ • \J 


0.003 


Jet Energy Resolution 


0.003 


0.003 


Other 


0.01 


0.010 


Total Detector 


0.02-0.05 


0.034 


Model uncertainty 






PYTHIA6-NLO 


0.00-0.05 


0.032 


Offset 


0.021 


0.021 


Scale 


+(0.01-0.05) 


+0.029 


-(0.01-0.02) 


-0.011 


PDF 


+ (0.002-0.004) 


+0.002 


-(0.002-0.007) 


-0.003 


MC Statistics 


0.005 


0.005 


Non-pert. Corr. 


0.002-0.014 


0.002 


Total Model 


+ (0.02-0.07) 
-(0.01-0.05) 


+0.044 
-0.034 


Total 


+ (0.03-0.09) 
-(0.03-0.08) 


+0.055 
-0.048 



we determine 95% CL limits on the contact interaction scale A. 

We summarize the determination of the limit in Fig. [3] We show TZll versus A for the data and 
for the SM expectation (with la and 2a bands) along with the highest value of TZll excluded 
at the 95% CL with the CL S method. The expected exclusion region comprises those values 
of A for which the SM-expected TZll (conditioned by the observed numbers of events n tot /) is 
less than the 95% CL S contour, and is seen to be A < 2.9 TeV. The observed exclusion region 
comprises values for which the measured TZll is less than the 95% CL S contour, and is seen 
to be A < 4.0 TeV. The observed limit is higher than expected because for m« > 1.4 TeV the 
measured R, 7 is lower than its expectation under the SM. 

In summary, we present a measurement of the dijet centrality ratio in 7 TeV proton-proton colli- 
sions. The dijet centrality ratio is found to exhibit little dependence on the dijet invariant mass 
and to agree with the expectation of the Standard Model. We exclude quark compositeness 
described by a contact interaction between left-handed quark fields at energy scales of A < 4.0 
TeV at the 95% CL. This is the most stringent limit to date. 
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Figure 3: Summary of the limit for the contact interaction scale A. We show 1Zll versus A for 
the data (solid line), the 95% CL S (dashed line), and the SM expectation (dotted line) with la 
(dark) and 2a (light) bands. 
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